Soliton molecules in trapped vector Nonlinear Schrodinger systems 
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We study a new class of vector solitons in trapped Nonlinear Schrodinger systems modelling 
the dynamics of coupled light beams in GRIN Kerr media and atomic mixtures in Bose-Einstein 
condensates. These solitons exist for different spatial dimensions, their existence is studied by means 
of a systematic mathematical technique and the analysis is made for inhomogeneous media. 
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Since the begining of its history, physics has studied 
simple objects and the way they arrange to form more 
complex structures. Some remarkable successes included 
the atomic theory of matter, the structure of nucleus in 
terms of protons and neutrons and the substructure of 
nucleons in terms of quarks to cite a few examples. 

Elementary robust objects made of light have been 
known since the 70's. In fact, spatial optical solitons — 
self-trapped states of hght with particle- like properties — 
have attracted a considerable attention during last years 
as possible building blocks of all-optical switching de- 
vices where light is used to guide and manipulate light 
itself ||. Another field where robust solitonic struc- 
tures have been recently found is that of Bose-Einstein 
condensation where the dilute quantum gas supports ro- 
bust structures such as one-dimensional dark solitons , 
bright solitons [Q or vortices 

In nonlinear optics, the robust nature of spatial opti- 
cal solitons allows to draw an analogy with atomic 
physics treating spatial solitons as "atoms of light" . Fur- 
thermore, when several light beams are combined to pro- 
duce a vector soliton, this process can be viewed as the 
formation of composite states or "molecules of light" . 

Several structures of this type have been studied pre- 
viously: dipole and multipole vector solitons |^, ||, ^ , 
self-trapped necklace ring beams ]lO{ , rotating propeller 
solitons ||ri| and rotating optical soliton clusters Q . 

In the field of Bose-Einstein condensation (BEG) vor- 
tex solitons arrange to generate lattice structures in ro- 
tating and non-rotating scalar systems In mul- 
ticomponent condensates simple stationary solutions in- 
volving some type of dynamical equilibria between the 
constituent solitons have been described |l^ . 

The purpose of this letter is to describe and ana- 
lyze in detail a way to build multidimensional "soliton 
molecules" in vector systems with applications to nonlin- 
ear optics and to Bose-Einstein condensed gases. As new 
features these solitons will be analyzed using a system- 
atic mathematical technique, the analysis will be made 
for inhomogeneous media and the soliton molecules exist 
for different dimensions D — 2,3 thus providing the first 
soliton molecules shown to exist for Z? = 3. 

The model.- We will consider a system of N complex 



fields ui{t, r),U2{t,r), ujv(^, r) ruled by the equations 



idt Uj{t, r) 



--A + V{r) + U,{t,r) 



Uj(t,r), (1) 



for j = 1, N. The coupling term is given by Uj{t, r) = 
J2k9]k\uk{t,r)\'^ with gjk G K. Eqs. (0) are a set of 
Nonlinear Schrodinger equations (NSE) which in BEG 
problems describe multicomponent systems, Uj being the 
wavefunctions for each of the atomic species involved 
[l^ , . In optics these equations describe the incoherent 
interaction between the slowly varying envelopes of the 
electric field in paraxial beams in Kerr media. We choose 
V{r) = r^/2 which corresponds to an isotropic magnetic 
trapping in BEG and to a GRIN fiber in the optical case. 
Single component case: Soliton atoms.- Let us first con- 
sider the scalar case {N = 1). Solitons or stationary 
solutions of Eq. in the scalar case have the form 
u{t,r) — (f)f^[r)e^^* and satisfy 
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For a fixed norm of the solution ||0||l2 :— J \(j)\'^dV, any 
solution to Eq. (^) will be a valid soliton atom for our 
purposes. The simplest case corresponds to a nodeless 
ground state solution which for strong interaction is close 
to the Thomas-Fermi solution (thus a quasi-compacton 
type of solution) and in the small interaction case (as 
it happens in nonlinear optics) it is close to a Gaussian 
function. Many other complex stationary solutions to 
Eq. (H) are possible such as vortices. These objects will 
play the role of "atoms" in what follows. 

The relevant property to be used here , is that any 
function of the form 

^.(r,t) = 0^(r-R(i))e'['^*+«("-^*)l, (3) 

is a solution of the scalar time-dependent NLSE provided 
R(t) satisfies ^R + R = and e{r,t) = (r, ^R)-|-/(t), 

where f{t) = (J^R)'^ — R^ dt. This means that ex- 
act scalar time-dependent solutions whose center evolves 
according to harmonic oscillator type equations and pre- 
serve the shape of the stationary solution during evolu- 
tion can be built. Remarkably, this property is not ex- 
clusive of the Kerr nonlinearity and it is valid for any 
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type of nonlinear term for which localized solutions ex- 
ist. However, if A'' > 2, the property described previously 
does not hold except for very trivial situations. Because 
of the cross-interaction. In a general situation, the pulses 
collide and lose their individuality. Here we want to give 
some ideas on how to overcome this problem and build 
stationary nontrivial vector solitons of Eqs. (|l|). 
Formalism for the multicomponent case.- Let us define 
the modulus, n^, and phase, (j>j, of each wavepacket 
through Uj = Y^exp {i^j). Let us define also the center 
of mass of species j, Rj(i) = / dV nj r and their total 
momenta Fj{t) — J dV rij'S/ipj, whose evolution laws are 



dt ^ 
d_ 

di 



(4) 

(5) 



The first term in the r.h.s of (|^) corresponds to the ex- 
ternal potential, while the nonlinear force is given by 



J dV {nkVrij 



njVrik) 



(6) 



Generally F^-fc may be a rather complex function, but if 
the Uj for j — 1, N, have small enough overlapping, it 
is possible to argue that these forces would be central. To 
prove this affirmation let us first notice that for the scalar 
case and far from the center of the wavepacket, the self- 
interaction is small and u can be described by the linear 
theory to be n (r) cx e^^ [r^™ + O (r^m-a^j f^j. ^ _^ qo. 
Similar considerations apply to the multicomponent case. 
If the wavepackets are separated, then in the region where 
the Uj overlap (i.e., the region which determines the value 
of the integral in (^)) the nonlinear terms will be small 
and 



Vuj = -2(r-Rj) 



from which it follows that 
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|r-R,| 



F,fe - 2 (R, - Rfc) J dV UjUk + O (^^^ 



Rfc 



(7) 



(8) 



i.e., if the wavepackets are separated, in the leading order 
the inter-mode force is central. 

As to the factor J dV nju^, let us evaluate it, for 
illustrative purposes, for the case when Uj,k are solu- 
tions of linear problem {gjk = 0) of the type Uj = 
TVj |r — Rj 1^™^ e"*-""^-^^^^ , where Nj are normalization 
constants (these functions include the fundamental mode 
rrij = and vortex states) 



dV rijUk = -NjNue~2\^^-^>'\ Kra„mA\^j - Rfcl 



(9) 

are polynomial factors, the lowest order 
ones for D = 2 being Xo,o(Q) = 1, Ko.i[Q) = 0V4 



where K, 
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FIG. 1: Potential W{Q) for gjk — (dashed line) and non- 
linear potential Wf{Q) (dotted line) for: (a) Interaction of 
two gaussian modes, (b) Interaction of a gaussian with a vor- 
tex mode. Solid lines represent the total effective potential 
for g/2TT = 20 (upper solid line) or g/2TT — —20 (lower solid 
line) (c) A situation with L = 0. Solid line: total potential for 
two gaussian modes with g/2TT = 5. Dashed line: potential for 
gaussian plus vortex modes with g/2TT = 20. The arrows mark 
the minima corresponding to non-rotating vector molecules 



and Ki^iiQ) = [2 ~ Q'^ + \Q^] /4. Thus, in the case of 
small nonlinearities, we can estimate the force between 
the wavepackets by using Eq. (|^). For our purposes the 
specific form of the interaction is less crucial than the 
fact that the forces Fjk are central. 

Two- component case.- For the two-component sym- 
metric {gi2 — 921 — g) case Eqs. (0) read 



(10a) 
(10b) 



where F = / dV (n2Vni — niVn2). The modified 'total' 
center of mass R = Ri-t-R2 and momentum P = P1-I-P2 
do not 'feel' the force F since they satisfy ^R = P and 



dt 



^P = — R, i.e. they oscillate harmonically. 

The most interesting quantity is Q = Ri — R2 (for the 
general multicomponent case we define Qij(i) — Ri(t) — 
Rj(i)), which gives the separation between the centers of 
mass of the two components and evolves according to 



df^ 



Q = gF 



(11) 
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FIG. 2: Evolution of initial data ui{r,t) = <^o(r- 
U2{r,t) = Mr + Ri(0))e-3"« ^j^j^ ^ ^ ^ 
J njdV = 1, Ri(0) — (3,0) and 0o is the scalar ground state, 
(a-c) Density plots of ni(r, t) + n2(r, t) on the spatial region 
[-8,8] X [-8,8]. (d) Evolution of Ri{t): Xi{t) (solid) and 
Yi{t) (dashed), (e) Evolution of Q(f). 



As discussed above, if Q is sufBciently large, the force 
F oc Q and it can be presented in potential form F = 
--^WfiQ) = -W{{Q)Cl/Q and then Eq. (|l|) reads 



(12) 



where Wtot — +gWi{Q). It is interesting to note 
that in the approximation of independent wavepackets 
the potential Wf is given by Wf{Q) = / dV nin2. 

It is evident that Eq. (|l^ has two constants of motion: 
energy E = ^(Q, Q) + W^tot(Q) and angular momentum, 
L = Q X Q. Since = + L^/Q'^, we can write E = 
iQ' + Weff (Q), where W,s(Q) = ^Q^ + J^^ + gWf{Q). 
There will exist an equilibrium distance Qo, for which 
the effective potential Wes is minimized. The reason is 
that the centrifugal contribution is divergent for Q ~ 
while the trap potential is unbounded for Q — > oo and 
the effective nonlinear interaction should decay for large 
Q and have a maximum finite amplitude [Fig. |l[ . 

In particular, if the nonlinear interaction is small then 
Qq" ~ L^/^ thus leading to a stationary rotating solu- 
tion of the vector system (||) provided the distances be- 
tween the components are kept large enough. The larger 
the nonlinear terms the larger will the deviation of Qo 
from Qo" be (Qeq > Qo for gjk > and Qeq < Qo for 
9jk < 0) [Fig |l|(a,b)]. When L = 0, the combination 
of the trap force and the nonlinear term may also have 
minima [Fig. |l|(c)] corresponding to non-rotating soliton 
molecules. However, due to the fact that the equilibrium 
distance is comparable to the wavepacket's widths, our 
result must be taken only as an indication that such soli- 
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FIG. 3: Evolution of a three-component system with 
uj{r,0) = (i/0f)e-((— .)'+(!'-!'j)')/2e'(-.j-+"«,!/)^ 
{X3,y3) = do(cos(27rj/3),sin(27rj/3)), {vxj,Vyj) = 
do(sin(27rj/3), — cos(27rj/3)) for j — 0,1,2. Parameter 
values: gij = 3,cio = 2, / rijdV — 1. (a-c) Density plots of 
X;^=inj(r,t) (d) Evolution of Ri: {Xi{t),Yi{t)) and Qi2{t). 
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FIG. 4: Same as Fig. |^ but with gtj — 5. (a-c) Density plots 
of X^'^^i '^j(i'i*)- (d) Long-time evolution of Qi2{t). 



ton molecules exist and not as a proof that they will be 
well approximated by combinations of scalar solitons. 
Examples of soliton molecules.- Let us now present sev- 
eral examples of the soliton molecules discussed previ- 
ously. First we have studied the case of a pair of weakly 
interacting soliton atoms. The results, obtained with a 
symplectic second order in time split-step integrator are 
summarized in Fig. |^ where it is seen how the small in- 
teraction induces only small harmonic oscillations of Q(t) 
without appreciable distortion of the wavepackets. 

In Fig. H three gaussian solitons interact more strongly 
due to the larger number of components and the smaller 
distance between the beams. In this case the oscillations 
of the distances between components, Qtj, remain small 
[Fig. 1(d)] ahhough some oscillation of the positions of 
the beams is appreciable [Fig. ^(a-c)]. Although this 
particular configuration is stable, the present example is 
a three-body problem for which many behaviors are pos- 
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FIG. 5: Same as Fig. ^ but with a gaussian soliton replaced by 
a vortex soliton. (a-c) Density plots of ^^(r, t). (d) Evo- 

lution of the inter-mode distances Qi2(t) (solid), and Qi3{t) 
(dashed) . 
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FIG. 6: 
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interaction of two gaussian solitons 

M2(r,0) = 



Ui(r,0) = (l/^3/4)g-({.-2)-+j;-=+.-=)/2g-2,y 

(l/7r3/4)e-«-+2)'+^'+-')/2e2»!' with g^^ = 10. (a-e) Isosur- 
face plots of ni{r,t) + n2{r,t) for ni + n2 = 0.01. (f) Evo- 
lution of the inter-mode distance Q{t) (upper solid line), and 
Xi{t) (lower solid line). 



sible: stable solutions, resonances, chaos, etc. In fact, if 
the values of the nonlinear coefficient are increased the 
beams deform and the inter-mode distances Q{t) suffer 
strong oscillations [Fig. ^(d)] although the structure re- 
mains stable after long periods of time containing about 
one hundred revolutions of the soliton around the center. 
Finally, in Fig. ||, it is shown that if one of the gaussian 
solitons is replaced by a vortex soliton the asymmetry of 
the interaction and the longer interaction range of the 
vortex soliton (for which n r^e"*" ) make this configu- 
ration unstable and the initial configuration is destroyed 
after a few rounds. 

We have also analyzed several soliton molecules in 
three spatial dimensions. In Fig. |^ we summarize the 
results for the situation of a stable configuration made of 
two weakly interacting gaussian solitons. The evolution 



for long times shows that the inter-component distance 
Q{t) suffers only small oscillations [Fig. ^(f)] which man- 
ifest on the plots where more interaction is apparent [Fig. 
|(c,d)]. 

Conclusions.- In this paper we have presented several 
soliton molecules built up from scalar solitons of the 
trapped NSE. The method presented here allows to gen- 
erate many different multidimensional soliton molecules. 

Between the possible applications of the vector solitons 
described here, in the field of optics they could be used as 
a way to transmit different laser beams in a single opti- 
cal fibers. In any case the vector solitons presented here 
represent another step into the comprehension of com- 
plex objects sustained by nonlinear forces with promising 
applications in different fields. 
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